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Abstract 

The neoclassical prescription to use an equation of motion to de- 
termine the electrostatic field within a tokamak plasma is fraught with 
difficulties. Herein, we examine two popular expressions for the equi- 
librium electrostatic field so determined and show that one fails to 
withstand a formal scrutiny thereof and the other makes a prediction 
for the electron temperature profile which does not compare well to 
that commonly seen in a tokamak discharge. Reconsideration of the 
justification for the presence of the equilibrium electrostatic field in- 
dicates that no field is needed for a neutral plasma when considering 
the net bound current defined as the curl of the magnetization. With 
any shift in the toroidal magnetic flux distribution, a dynamic electric 
field is generated with both radial and poloidal components. 



1 Introduction 



The neoclassical prescription to use an equation of motion to determine the 
equilibrium electrostatic field within a tokamak plasma, rather than Gauss's 
law or Poisson's equation [HOE], is fraught with difficulties alleviated by 
insisting on the neutral fluid limit rather than the quasineutral approxima- 
tion [H [5] . Herein, we examine two popular expressions for the electrostatic 
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field so determined, the first established from electron momentum conserva- 
tion and the second from the Pfirsch-Schliiter current, and show that one 
fails to withstand a formal scrutiny thereof and the other makes a predic- 
tion for the electron temperature profile which does not compare well to that 
commonly seen in a tokamak discharge. Reconsideration of the justification 
for the presence of the equilibrium electrostatic field indicates that no field is 
needed when all bound currents are considered and that the field must either 
vanish or result from sources for a model based upon the physics of a neutral, 
conducting fluid. Shifting of the toroidal magnetic flux density results in a 
circulating electric field with both radial and poloidal components. 

At issue is the validity of the quasineutral approximation, which allows for 
a divergenceful electric field in the absence of a non- vanishing space charge 
density p e = J2 S n s e s, formally expressed as V ■ E ^ for p e = 0, and 
requires the determination of the electrostatic potential from an equation 
of motion. However, the quasineutral approximation does not respect the 
mathematics of electrodynamic field theory, 

V • E = Pe /e . (1) 

From a particle physicist's field-theoretic perspective [SI 13 El EJ [10] , the gauge 
invariant Maxwell field tensor F^ v = d^A v — d v A 11 is known to carry only 3 
scalar degrees of freedom in media, not 3 for each of the electric and mag- 
netic fields, embodied by the four-potential = ($/c,A) subject to the 
gauge condition and coupled to sources given by the conserved four-current 
J M = (cp e , J) through the inhomogeneous Maxwell equations d u F^ u = poJ^, 
and the homogeneous equations are recognized as the Bianchi identity for 
electromagnetism given by the field equation for the dual tensor d v F^ u = 
and are satisfied identically when written in terms of the electromagnetic po- 
tential hence do not determine any degrees of freedom, thus the electrostatic 
field is determined by the space charge density p e and not by an equation 
of motion. The equations by Maxwell may be expressed succinctly using 
intrinsic, geometric notation as d *d A = J in terms of the exterior derivative 
d, the Hodge dual *, the connection 1-form A, and the current 3- form J, 
as given in many standard quantum field theory texts, such as Ryder [S], 
or more esoteric monographs, such as Davis [7]. What this shows is that 
Gauss's law may not be isolated from the remainder of the source bearing 
Maxwell field equations, V x B — p e dE/dt = PqJ, as they are but one 
unit of truth. 
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Note that we are not criticising the kinetic approach to plasma calcula- 
tions, which in its original inception as the Vlasov- Maxwell system [TT] fully 
respects the microscopic electrodynamic field theory, but rather the neoclas- 
sical (non-classical) fluid model based on the quasineutral approximation, 
which does not respect the macroscopic electrodynamic field theory. While 
various approximations are made in the drift kinetic equations [T21 [T5] , most 
practical numerical evaluations address Poisson's equation directly [Hj. With 
this manuscript, we examine in detail the mathematical difficulties one en- 
counters when following the neoclassical prescription. Similar discussions on 
the existence of the whistler oscilliton in geophysical plasmas [151 EE E3 EH] 
and of the helicon wave in propulsion devices |19 : , ;2Dj are noted. Criticism of 
the quasineutral approach in a cosmological context [2TJ [22] has also recently 
appeared, and it is time for the fusion community to address these difficulties 
head on. 

2 Electrostatic field from the species equa- 
tions of motion 

Some geometric nomenclature and vector identities are given in the Ap- 
pendix. The first equilibrium d/dt — > electrostatic field under consideration 
is one commonly used in the analysis of tokamak experiments [23J [2J], E5J EH] , 
determined by integration of the electron poloidal and ion radial equations 
of motion in the large aspect ratio, concentric circular flux surface approx- 
imation. The use of concentric circular flux surfaces with a toroidal inte- 
grating measure is pursued herein to remain consistent with the model as 
presented in the literature, as are the expansions of the electron density 
n e = n° e {r) [1 + n c e (r) cos 6 + n s e {r) sin 6] and electrostatic potential. Refer- 
ence [23] states that "the electron momentum balance can be solved for 
(pc,s = (£c,s j £ _ n c,s j £ (^ e ^fi /T e ), which represents the poloidal asymmetry in 
the electrostatic potential." Let us examine that statement in detail. 

This model [23], EU E5] writes the equilibrium poloidal equation of motion 
for arbitrary species s as 

[n s m s (V 8 ■ V) V s + V ■ II S ] ■ 9 + dp s /rd6 - F s9 + n s e s {V ar B^ - E e ) = , 

(2) 

where p s = n s T s for T s <— ksT s and F s is the friction term, and takes the 
poloidal component of the electrostatic field on a flux surface at r in Coulomb 
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gauge as 



d$(r,6) 
rdO 







rdO 



<P°(r) 



i [1 +<P c (r)cos9 + $ s (r)sm9] (3) 
[<P S (r) cos 9 -<P c (r) sin 9] , (4) 



where $ is the electrostatic potential, indicating an expansion around <£°(r) = 
— f drE® ^ for a last closed flux surface at r = a, where the radial 
electrostatic field is calculated from an ion equation of motion [22]. The 
resulting evaluation of the flux surface unity, cosine, and sine moments of the 
electron poloidal equation of motion with dT e /d6 = (where other terms are 
assumed negligible at equilibrium), 



T e dn e /d9 = —en e rE g , 



(5) 



defined by the expressions (A)^j,c,s} = ^e£0{l,cos0,sin0}(l + e cos9)A/2ti, 
yields three equations which have only trivial solution. Specifically, we have 
the system of equations 



U 
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S 



en s e T e 
n c T P 



e<P°(e& s + n c e & s - n s e $ c ) , 
e0°(40 s + 3enl<P s - en s e $ c )/A 
e$°(4<P c + en c e <P c - en s e <P s ) /4 , 



(6) 
(7) 
(8) 



valid V e, n% n s e , and T e . The terms with factors of e = r/R above are strictly 
due to the toroidal geometry and would disappear for a cylindrical plasma 
column R — >• oo; one cannot address the extension to toroidal geometry 
of other aspects of the model [27] without addressing the extension here. 
Solution in pairs given finite (fixed) <P° yields inconsistent values of (£> c ' s and 
an overdetermined system, which therefor has no solution, thus the poloidal 
electrostatic field in this neoclassical model, which fails to consider the 0(e) 
terms within the cosine and sine moment equations, is unphysical. Failing to 
include the 0(e) terms indicates expressions applicable only on the magnetic 
axis r = 0, where e = for Rq ^ oo, yet this neoclassical model is commonly 
used to address the physics near the edge of the confinement region [23J. 
This system may be put into linear, homogeneous form Ax = by dividing 
through by 
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thus its only exact solution for e ^ is trivial (1 / <P° , <P C , <P S ) = (0,0,0), 
which we interpret to mean exactly what it says, that it is solved when 
<?°(r) = ±oo, displaying its unphysical definition when determined from 
the electron poloidal equation of motion. The matrix has rank 3, thus the 
solution has free parameters [28]. For a cylindrical column e — > 0, one 
recovers a matrix of rank 2 and the solution $ c,s = nl ,s (T e / e<P°) with <£>° a 
free parameter. If one assumes that the smallness of the coefficients may 
be represented by e, eg n c f = en^, then the offending terms near the edge 
where e approaches 1/2 represent up to a 19% correction to the leading order 
equations. As non-vanishing <P C,S are an integral part of the development of 
this neoclassical model and appear in its remaining equations without a pref- 
actor of <P° through substitution, the validity of its conclusions is in jeopardy. 
Note that a putative non- vanishing radial electrostatic field without poloidal 
variation demands the existence of no poloidal electrostatic field, else the 
poloidal variation to the potential ruins the poloidal symmetry of the radial 
field; if that radial field is determined from a radial equation of motion then 
the associated poloidal field is determined by the poloidal dependence of that 
equation. 

One might think to alleviate the difficulty by invoking the logarithmic 
derivative, writing the electron poloidal equation of motion as <9(lnn e — 
e@/T e )/d9 = 0, with solution C(r) = n e (r,6)exp[-e&(r,6)/T e (r)]. Expand- 
ing the exponential gives, for — oo < e<P/T e < oo, 

e<P 1 fe<P\ 2 1 (e&\ 3 

(10) 

and to each order in e@/T e , taking the flux surface moments yields three 
equations to solve. Using the expansions for n e and $ above, the zeroth or- 
der equations give C{r) = n° e and rff = 0, and the first order equations yield 
e-dependent solutions which, upon restriction to the magnetic axis e — 0, are 
C{r) -> n°[l - K) 2 /2 - (^) 2 /2](l - e<P°/T e ) and <P C > S -> n c e > s (T e /e$° - 1). 
Agreement is found with the previous <& c,s = nl' s (T e /e^>°) provided e<P°/T e 
1. Considering now a ratio e<P°/T e ~ 10 ^> 1, such as may be found in the 
plasma core [2U [29J, [30], the expansion in Equation fflOl) requires 10 orders 
before the absolute magnitude of the terms starts decreasing, indicating that 
a low-order approximation is not an exact solution. Considering the second 
order equations leads to solutions for C(r) and <P C,S only when e = which 
diverge when n^ s = 0, and no solution has yet been found to the nonlin- 
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ear equations for general e. The difficulty with <P° encountered above has 
simply been shifted to the "undetermined" function C (r) , which is perfectly 
determinable in principle from the system of equations. We remark that, 
assuming the model's system to be well-determined (ie equal numbers of in- 
dependent equations and degrees of freedom) before invoking the logarithmic 
derivative, the introduction of C (r) without an additional equation for its de- 
termination leaves the system under-determined — if C(r) is identified as an 
independent degree of freedom, then it needs its own independent equation, 
else if not, then one of those degrees of freedom must be determined by the 
unity moment equation here. Pursuing the argument, as the electron density 
71° divides out of Equation (j5J), its poloidal variations n c f are determined by 
continuity {dn e /dt + V • n e V e — h)c,s — for particle source rate n, and 
the thermal energy T e is determined by the heat equation and given as in- 
put from experimental measurement for the analysis, the remaining degree 
of freedom <P° must be determined by C(r) as given by the poloidal force 
balance system, Equation (jUJ). 

Notwithstanding the above, let us now consider an alternate expansion of 
the exponential which extracts the dependence using £ = <P C cos8 + <P s sin 0, 
so that 

C(r) = n e exp(- e ^ /T e )(expO" e< ' 0/T ' ! , (11) 
« n e exp(-e<Z>°/T e ) (l + £ + £ 2 /2 + . . .) _e * ° /Te , (12) 
« n e exp{-e<P°/T e ){l + cos + n* sin 0) x (13) 
e<P° 1 fe<P°\ 2 2 1 

c+aCsr) e* + -J ■ (14) 

noting the use of low-order truncations to the two infinite series (which can 
never be as exact as the algebraic method above) and that technically the 
expansion found in Equation fflOl) is more democratic in its treatment of 
the powers of e<P°/T e . Taking the flux surface moments as before, we find 
nearly the same first order solutions, whereupon reduction to the magnetic 
axis (or cylindrical geometry) yields C(r) — > n°exp(— e$°/T e )[l — {n c e ) 2 /2 — 
(rig) 2 /2] and <P C,S — > nl' s (T e /e<P°), and for the second order equations the 
solution on the magnetic axis diverges and for general e cannot be found. 
Inserting factors of e at the cost of the loss of the cylindrical interpretation 
of e — > 0, which still identifies the magnetic axis in toroidal geometry, yields 
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Equation fllip in terms of the tilded coefficients (which need not be small), 
C(r) w n° e exp(-e0°/T e ) [1 + e(n c e cos 9 + n s e sin 9)} x (15) 



—e($ c cos # + $ s sin #) + i — e 2 (<Z> c cos + <Z> S sin ^- 

J- e. \ e / 



,(16) 



and affords additional control over the terms retained in the expansion. The 
resulting system is of the form = Uo + U 2 e 2 + U4E 4: = Ci£+C 3 e 3 = S1E+S3E 3 
and has not solution until reduced to order 0(e 2 ) which, as a factor of e 
cancels out of the (. . .)c,s moments, is the first appearance of an explicit e 
dependence in the equations. Its expression for <E> C contains terms with l/e 2 , 
thus the physical <P C —>■ 00 when e — 0. Only when each equation is reduced to 
its leading term does the system return the solutions C(r) = n° e exp(— e<P°/T e ) 
and <P C,S = nl' s (T e /e < P°). Comparing the leading coefficients in Equation <TT~0T) 
gives the ratios 1:10s: 55e 2 for e<P°/T e ~ 10, and allowing for &P° /T e ~ 1 
yields 1 : e : e 2 which still represents up to a 25% correction. The lesson here 
is that knowledge of a series' convergence need not imply that a low-order 
approximation is an adequate representation of the function. Concern over 
the expansion in no way detracts from the observation that the physics of 
the situation is embodied by the algebraic Equations ([6]) which contain no 
exponential factor thus represent a more exact method of solution. 

Finally, we consider the consistency of using two different species' equa- 
tions of motion to determine the species independent electrostatic potential, 
taking the radial electrostatic field from the equation of motion for arbitrary 
ion species j as 

1 Op ■ 

E r = 7; 1- VsjBo — VgjB^ , (17) 

HjCj or 



and dropping the convective term as is standard practice in the field [24 
[26] . The usual evaluation of that expression from experimental measurements 
neglects any poloidal dependence; however, when the intrinsic variation of 
quantities induced by the geometry is considered [3T], given in the large 
aspect ratio e < 1, concentric circular R r = R , flux surface approximation 
B = (0, B e , B$) by 



B = B°/(l+ecos9) , V 0j = V e ° J /(l+ecos9) , = Vj^l+ecosfl) , n,- 

(18) 

where A is the average of the values of A on the vertical midplane, a geo- 
metric dependence is introduced. Using these values (r dependence implied), 
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we find 



E r (r,9) 




+ V« ] B*-V* J Bl/{l + ecos6)\ 



(19) 



for species pressure gradient p'- = drijTj/dr, thus E r {6) ^ E®. 
the denominator reveals a power series in ecos9, 



Expanding 



E r (r,9) 




(r) + E;(r) cos 9 




+ - V^Sj (1 - 2e cos 9 + 3e 2 cos 2 - 4e 3 cos 3 




(22) 



(21) 



and as integration with respect to r does not affect the 9 dependency, the 
potential associated with the radial electrostatic field relative to its central 
value, $£ r (r, 9) = — J Q r drE r (r, 9), may be written as a cosine series, 



where $ c Er = drE}/ j^drE® ^ in general, noting that the potential on 
the last closed flux surface at r = a is not single valued. Returning now 
to the electrostatic potential appearing in the poloidal equation of motion, 
expressed to leading order as 



where <E> E/j (r) = — J a drE®. ^ is the potential relative to that of the last 
closed flux surface, with solution <E> c Ee = n^ ,s {T e /e^>). No loss of generality 
ensues if one redefines the potential relative to its central value. With vanish- 
ing extrinsic poloidal dependence to the electron density, n^ s —>■ as above, 
the potential retains no explicit poloidal dependence, $E e (r, #) — > <E> Efj (r). 
Thus, we conclude that the electrostatic potentials associated with the ra- 
dial and poloidal electrostatic fields evaluated from the ion and electron equa- 
tions of motion by this neoclassical model are inconsistent, $E r 7^ ^E e , as 
Equation (!23j) does not equal Equation (1241) in the case of vanishing density 
asymmetries * = 0. 



$ Er (r, 9) « $° Er + <P Er cos 9 = <P° Er (1 + $ Er cos 



0) 



(23) 



$ Es (r, 9) = $ Ee (r) [1 + <P Eg (r ) cos 9 + (r) sin 9} , 



(24) 
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3 Equilibrium field from the Ohm's law equa- 
tion 



3.1 Derivation 



Examining the expression of another leading contender for the equilibrium 
electrostatic field [3] evaluated from the Ohm's law equation and the Pfirsch- 
Schliiter current, one may put its poloidal component 



En 



{E ( j ) B ( j > /B g )B 2 E^B 4 



(B 2 /B d )B e B 6 



, ( (l/B e )B 2 



1 

B~ f 



(25) 



into the form 

E e = Eg [2e cos 9 - (e 2 /2) cos 29] , (26) 

when the Shafranov shift is neglected, as in the concentric circular flux surface 
approximation of above, upon application of Stokes' theorem to Faraday's 
law [32J, ie by requiring § d9Eg = 0. Note that this neoclassical model for the 
poloidal electrostatic field differs distinctly from that of the previous section 
in detailed functional form. Inserting Equation (1261) into Equation (jSJ) and 
taking the flux surface Fourier moments yields three equations which have a 
nontrivial solution only when expanded to order 0(e 3 ), given by 



" n c e ' 




K 









e 3 /(6e< 



±eV3e 4 - 168e 2 + 192/(18£ 2 - 24) 
±4(T e /e J R )Av / 3£ 4 - 168e 2 + 192 



(27) 



thus the presence of a poloidal electrostatic field of that form should be 
accompanied by a potentially measurable shift in the electron density profile. 
Note that the derivation immediately preceding is slightly inconsistent, as the 
associated electrostatic potential <P axi = <P (2esm9 — (e 2 /4) sin 29) is of the 
correct harmonic form [331 134"1 1351 [36] for axial geometry, as is easily verified 
in (Z = —r sin 9,R = Rq + t cos 9, z) coordinates via application of the axial 
Laplacian V 2 „ = d 2 /dZ 2 + d 2 /dR 2 to <P axt = <P Z[(R - R )/2R 2 - 2/R ], 
yet the flux surface average is done in toroidal geometry. Note that this <Po 
is not the of the preceding section but is a unit bearing constant which 
sets the scale. 

In order to achieve the correct harmonic form for tokamak geometry, 
the potential must satisfy the toroidal Laplacian, which in cylindrical co- 
ordinates E3 EH] is given by V 2 or = V 2 axi + d/RdR (note that the 
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expression for A* given by Hopcraft [2] is not the toroidal Laplacian and 
differs by the sign of the additional term), and the form of the geometric 
term hints at the solution. Direct integration yields the toroidal potential 
$tor = $axi(R ->• ln-R), from which E z = -$ [(\nR- R )/2R 2 - 2/Rq] and 
Er = —<PqZ/2RRI, noting that the introduction of the logarithm breaks the 
usually obvious relation between the symbol for the magnitude of a quantity 
and the units associated with that quantity — carefully pulling the units be- 
side the leading coefficients of expressions ensures that they are respected. 
From these, we determine the poloidal field to be 



~ r 06 \ r J \d0 dZ dO dRJ &o 6 



$o $o 



(28) 

where we identify E^ e = —<Po/r, and the corresponding radial field is 



Z^-(R-Ro)^' 

00 ^0 



= _d<P _ f-$o\ r fdZ d dR d \ <P _ 
dr \ r J \ dr dZ dr OR J <P 

(29) 

In (r,9,<p) coordinates, we have <P ior = ^o^sin^f— \ia(R + r cos 9) + 5Ro]/2Rl, 
and 



E ^-Sk\ r sme[ln ^ +r ™ e) -^ ]+r ^^'> ■ (31) 

As this is an electrostatic field within a neutral medium (ie one for which 
the net charge on a differential volume element vanishes) at equilibrium, 
Maxwell's equations V x E = and V ■ E = are satisfied. 

As electrostatic fields necessarily require a supporting charge density [6j 
S2J [391 El E], one may very well ask where these charges are. As the solution 
to Laplace's equation is uniquely determined by the boundary condition, we 
will find them on the boundary of the region under consideration, which in 
this case is the R/Rq weighted circle representing our outermost flux sur- 
face at normalized minor radius r/a — 1 upon collapse of the toroidal di- 
mension, giving us an inverse Dirichlet problem, which is usually defined as 
solving for the potential given the boundary condition. The source of these 
charges is intended to be the accumulation resulting from the divergence of 
the pressure gradient driven diamagnetic current in toroidal geometry, and 
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with —dp e / dt = V • J ps oc V • E, we find the rate of surface charge accumu- 
lation & b oc —f ■ E{r = a), noting that quantity represents the charge density 
on the weighted circle, so to get the physical charge density on a poloidal cir- 
cle of the toroidal flux surface we take &b — > dbRo/R. The Pfirsch-Schliiter 
current is supposed to flow along the electrostatic field at equilibrium to 
cancel this charge accumulation. 

From the expression for the poloidal field, one may determine a prediction 
for the associated electron temperature profile by solving the flux surface 
Fourier moments of the poloidal equation of motion as in the previous section. 



With electron density n e 



'l+n c e cos 9 + n s e sin 9) and temperature T e , we 







K 




. E e_ 





solve the unity, cosine, and sine moments (T e dn e /rd9 + n e eE e } UjC ,s — for 
(the pure numbers) n c e1 n s e1 and Eg, finding 

3e 3 /G 2 

±y/EeG 1 /2{F + 2)G 2 , (32) 
±8VQT e /e£G 1 

where F = 5Rq — \nR Q , in terms of the functions G*2( r ) = [{6F — l)e 2 — 8F] 
and 

Gi(r) = ^(F + 2)[{16F[4F - (1 + 3F)e 2 ]} + 9e 4 ] , (33) 

from which we obtain <& = ^8\/6R T e /eGi. Note that the quantities ap- 
pearing in, between, and following Equations (1321) and (133|) carry no units, 
as the units cancel out of the system of equations solved. The factor of Rq 
appearing as a scale factor in <P has a physical interpretation — for an equiva- 
lent supporting charge density and minor radius, as one scales the system by 
Ro, the magnitude of the charge distribution on the weighted circle increases 
as the major circumference of the torus, thus increasing the magnitude of 
the potential by a numerical factor of Rq. As <Pq is constant V r, the ex- 
plicit r = eRo dependence of Gi must cancel the implicit r dependence of 
T e (r), allowing us to write T e (r) /Gi{r) = T e (0)/Gi(0), thus Gi determines 
the r dependence of the electron temperature as a prediction of the hypoth- 
esis that the electrostatic field in a tokamak is of a form similar to that of 
Equation (1261) . Note that a similar compatibility relationship obtains from 
the bottommost of Equations (12 7\i given by the radical factor. 



3.2 Evaluation 

The following presents the results of a numerical evaluation of the expressions 
in the preceding subsection. The scale is set by Rq = 2 and T e (0)/e = 1, and 
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Figure 1: The electrostatic potential $ tor for toroidal geometry in arbitrary 
units seen from multiple views (a)-(d). The magnitude of the associated 
electrostatic field in arbitrary units seen from multiple views (e)-(h). 



r /a is the normalized minor radius. Some useful geometry is presented in the 
Appendix. The toroidal electrostatic potential <P tor is shown in Figure H](a)- 
(d) — note the correct harmonic form of a stretched circular membrane (eg 
a drum-head put on a hoop given by the boundary condition), consistent 
with the physics of potential theory [531 El] , and the obvious lack of bumps 
or poles, consistent with a neutral medium (the term "pole" actually means 
something physical [32j [3H [35], namely the effect on the potential of an 
isolated flux-source which ultimately is quantized in units of the charge). 
We stress that these are multiple views of the same single object and that 
the boundary is determined by the potential and not vice versa. In Fig- 
ure [Tfe)-(h) we display the magnitude of the associated electrostatic field, 
and in Figure |2](a)-(d) we display the electrostatic field and the associated 
supporting charge density. 

The normalized electron temperature profile in a tokamak often compares 
favorably (from unpublished results of the analysis behind Reference [24J) 
with the formula T e (r)/T e (0) ~ (1 — r 4 ) 2 , which misses the effect of the 
pedestal Tf ed /T e (0) ~ .2 but works surprisingly well for r/a < .8, and so that 
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Figure 2: The electrostatic field and associated supporting charge density 
in arbitrary units seen from multiple views (a)-(d). Comparison (e) of the 
neoclassical prediction for the normalized electron temperature profile, G\, 
versus an approximation to that found in tokamak experiments, (1 — r 4 ) 2 . 



is the assumed profile to which we will make our comparison in Figure EJ^e). 
The normalized predicted profile for the temperature is Gi(r)/Gi(0) and 
fails to fall sufficiently quickly to match the shape of the expected temper- 
ature profile. Thus, while comparison with our approximate formula hardly 
qualifies as a true experimental test of the hypothesis, the existence of an 
electrostatic field of form similar to that of Equation f[2T)j) is not supported 
by this analysis, even allowing for a pedestal at 20%. The observed electron 
temperature at the edge of the confinement region would have to exceed 90% 
of its central value to be compatible with the model of this section. 

Regardless of the preceding argument, relying as it does on an abbrevi- 
ated form of the electron poloidal equation of motion, the supporting charge 
density above is reminiscent of the surface charge on a current carrying con- 
ductor, a topic of recent interest [HH HH H21 US Hi]- Note that these authors 
rely extensively on Poisson's equation, ie Gauss's law. Despite the poloidal 
appearance of the initial expression, Equation fl25|) . the resultant electric field 
is primarily vertical and should lead to an E x B drift force helping to bal- 
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ance the magnetic pressure force in the major radial direction. Investigations 
are proceeding along these lines to determine whether the return flux of the 
fringing field from the finite central solenoid [15] can induce such effect. 

The motivation for the existence of the Pfirsch-Schliiter current leading 
to the poloidal electric field of this section is the cancellation of charge ac- 
cumulation arising from the pressure gradient driven diamagnetic current 
Jv P = — Vp x B I B 2 in toroidal geometry [3J, V tor • Jy p = —dp^ p /dt ^ 
for 

V- Jvp ^V.(^) (34, 
= Vp.(VxB)- B .(VxV rt +(BxVp) , v _^ (3s) 

where the final term is nonzero due to the poloidal dependence of B. How- 
ever, Jvp is but one component of the total diamagnetic current, properly de- 
fined as the curl of the magnetization J ^ a = V x M where M = —(p/ B 2 )B, 
which includes the effects of the pressure gradient driven current as well as the 
curvature and V5 drift currents jl6] and remains divergence-free regardless 
of the geometry V • J ma = 0, thus there is no space charge accumulation and 
no motivation for a cancelling current. (The common procedure of adding 
the particle drift currents to the fluid diamagnetic current we feel represents 
an over-counting of the underlying phenomenon, as all currents are the curl 
of either an H or an M.) The error here lies in not fully distinguishing the 
free and bound charges and currents as they appear in the Heaviside notation 
of the Maxwell equations and in trying to model a fully ionized medium as 
both a conductor and dielectric at zero frequency. 

4 Dynamic electric field 

Relaxing our requirement of equilibrium slightly, we consider now the effect 
of a shift in the magnetic flux density when the total magnetic flux remains 
constant. As Gauss's law is inviolate, any electric fields within the neutral, 
conducting medium of a tokamak plasma necessarily are driven by changes 
in the magnetic flux density at some location in space, giving us 

V • E = , V xE = -dB/dt . (36) 

These changes may result for fixed by a shift in the pressure distribution 
leading to dM^/dt ^ 0, as B//i = (H — M)h. Noting the similarity to 
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1.5 2 2.5 -0.5 0.5 

R Z 

Figure 3: Electric field generated from changing magnetization for param- 
eters given in the text, (a) Contours of the flux function follow the elec- 
tric field, (b) Magnitude of the electric field along the vertical midplane at 
#0 = 2. 



the laws of Thomson and Ampere allows one to define an electric vector 
potential E = V x F in the analogue of Coulomb gauge V • F = 0. Then 
V 2 F = dB/dt has the solution [25l H7J HH 09] for a circular loop source at 
(Z ,R ) given by 



K(k)-E(k) 



(37) 



expressed in terms of the complete elliptic integrals [50] with parameter k 2 = 
2b/ {a + b), where a = (Z - Z ) 2 + R 2 + R 2 and b = 2RR and its magnitude 
F° = —(dB ( f ) /dt)A 2 R /4:7i is supposed constant over a differential area A 2 . 

On a uniform grid (Z,R) in meters with spacing A = .01 m, we wish 
to evaluate for sources with opposite polarity around (±Z , R ) for Z = 
0.25 aligned to the vertical midplane at i?o = 2. Expressing the plasma 
magnetization [51] as M/H = (1 — -y/l — 4p/ fi H 2 )/2 for h = 0cosC + <9sin£ 
and supposing the pitch angle ( = eir/2 for e = Z /Rq and dH/dt = allows 
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one to write 

dB^/dt = -(cos()(H 2 - 4p/i2 )- 1/2 dp/dt . (38) 

Using some typical parameters n = 10 19 m -3 , T e = 3 keV, T, = 9 keV, 
B<f, = 2 T, B e = 0.2 T, and supposing the rate of change to the pressure 
is 1% per millisecond, lets one give a numerical estimate to F® hence the 
magnitude of the electric field. From Gauss's law one defines the flux funtion 
X such that 

* x = 4> x Vx = RE = RV x F , (39) 

whence \ — ~RF<f> upto an unphysical constant, where >3& is the contour op- 
erator |33]. The results of this evaluation are shown in Figure [31 where one 
sees that an electric field with magnitude approaching 1 mV/m for these pa- 
rameters may arise on the vertical midplane. Pursuing the analogy between 
E ~ B and J ~ —dM/dt, a force of repulsion 

F + _ = -^t Q {dM ± /dt) x E T (40) 

should appear between line sources of opposite polarity by action of the 
macroscopic Lorentz force [SU \52\ where the subscript indicates the 
source location, thus a perturbation to the pressure does not necessarily 
coalesce. 

5 Conclusions and outlook 

From the preceding analysis, we find that the use of an equation of motion 
to determine the equilibrium electrostatic field leads to inconsistencies, ei- 
ther internally or with expected experimental measurement of the electron 
temperature profile. Part of the problem lies in treating the poloidal and ra- 
dial components of the field separately, when there is only one electrostatic 
potential from which both components may be determined. The remainder 
lies in the neglect of Gauss's law popularly established within the plasma 
physics community rjlElEliaElEZlIMlESlESl^E^En], which relegates 
the defining relation for the electrostatic field to a position of subsidiary mo- 
ment. Any model which treats the plasma as a neutral, conducting fluid, 
where neutrality is understood to hold down to some scale smaller than the 
differential volume element used to define the continuum quantities, needs to 
respect all of Maxwell's laws, which are manifestly Lorentz covariant. Taking 
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a field theoretic perspective implies that the electrostatic field within a neu- 
tral medium must either vanish or result from sources and be supported by 
dielectric polarization. In conclusion, we determine that we are unsatisfied 
by either neoclassical model for the equilibrium electrostatic field within a 
tokamak plasma on the grounds that they do not respect the mathematics 
of field theory. An alternative description in terms of shifting magnetic flux 
may possibly account for any direct experimental observations of an electric 
field within the device [60, 61J. The thrust of our argument is simply that, 
if there is an electric field within a tokamak, then it must play by the same 
rules as every other electric field found in nature. 



A Appendix 

There are (at least) three useful sets of coordinate axes to describe a toroidal 
magnetic confinement device with concentric circular flux surfaces, namely 
(Z,R,(j>), (r, 6, <f>), and (f, J_,||), and in the infinite aspect ratio limit 
(R Q — > oo) we have the axial coordinate axes (Z,R,z) and (f,6,z). Note 
that the term "toroidal coordinates" means something very different to a 
mathematician than those commonly applied to a tokamak, which we call 
"tokamak coordinates." For a plasma with coaxial applied electric and 
magnetic fields and free current driving a circulating field, we note that 
(f,6,z) = (-ExB/EB,-Ex(ExB)/E 2 B,E/E) and (f,i, j|) = (-Ex 
B/EB, —B x(Ex B)/EB 2 ,B/B). Cylindrical coordinate labels (Z,R,<f>) 
relate to tokamak coordinates (r, 9, 0) via Z = —r sin 9 and R = Ro + r cos 9 
in the concentric circular approximation, where f, 6, and (p gi ye the radial, 
poloidal, and toroidal directions, respectively. The outermost minor radius 
of the confined plasma, given in meters by a, defines the normalized minor 
radius r/a, and R a is its centroid. The magnetic field B and current density 
J lie in isobaric surfaces given by Vp = J x B for a stationary equilibrium, 
defining the "flux surface" at radius r. In general, the nested flux surfaces 
are neither circular nor concentric. The relationship between vector compo- 
nents in the tokamak coordinates (Z, R, </>) <— (f, 8, </>) <— (f, J_, ||) may be 
succinctly expressed by 



" F z ' 




Fr 









— sin 9 — cos 9 
cos 9 — sin 9 
1 



1 
b e 
—bo b rh 





' F r ' 




F ± 




F \\ 



(41) 
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where || = B/B = (0, 60,6^). Various useful relationships are 



Z = — rsin# , dZ/dr = Z/r , 
R-R = rcos6 , dR/dr = (R- R )/r 



dZ/dO = — 
dR/dO = Z 



(R — R ) , 



(42) 



and for the logarithm, we have d In R/dR = 1/R where the units on the left 
are carried by the differential operator and the units on the right are carried 
by the result, which shows that the In R is a pure number which carries no 
units. 
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